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The half-filled Hubbard model on the Bethe lattice with coordina- 
tion number z = 3 is studied using the density-matrix renormaliza- 
tion group (DMRG) method. Ground-state properties such as the en- 
ergy E, average local magnetization (S z ), its fluctuations (S 2 ) — (S z ) 2 
and various spin correlation functions {S z (i)S z (j)) — (S z (i)){S z (j)) are 
determined as a function of the Coulomb interaction strength U/t. 
The calculated local magnetic moments (S z (i)) increase monotonically 
with increasing Coulomb repulsion U/t forming an antiferromagnetic 
spin-density-wave state which matches the two sublattices of the bi- 
partite Bethe lattice. At large U/t, (S z (i)) is strongly reduced with 
respect to the saturation value 1/2 due to exchange fluctuations be- 
tween nearest neighbors (NN) spins (\{S z (i))\ ~ 0.35 for U/t — * +oo). 
(S z {i) 2 } - {S z (i)) 2 shows a maximum for U/t = 2.4-2.9 which results 
from the interplay between the usual increase of (S z (i) 2 } with increas- 
ing U/t and the formation of important permanent moments (S z (i)) 
at large U/t. NN sites show antiferromagnetic spin correlations which 
increase with increasing Coulomb repulsion. In contrast next NN sites 
are very weakly correlated over the whole range of U/t. The accuracy 
of the DMRG results is discussed by comparison with tight-binding ex- 
act results, independent DMRG calculations for the Heisenberg model 
and simple first-order perturbation estimates. 



I. INTRODUCTION 

Bethe lattices or Cayley trees have often been the basis of very attractive models 
for the theoretical study of various properties of solids. A Bethe lattice is completely 
characterized by the number of nearest neighbors z and by the lack of closed loops. 
The later feature simplifies calculations considerably allowing in many cases to ob- 
tain useful insights on the physics of complex problems, for example, in theory of 
many electron systems or in the theory of alloys and other disordered systems. Re- 
cently, the interest in the study of strongly correlated fermions on Bethe lattices has 
been renewed by the advances achieved in the limit of infinite spatial dimensionality 
(d — > oo)EJ. Indeed, this lattice provides a systematic mean of realizing the d = oo 
limit by letting z — >■ oo and scaling the nearest neighbor (NN) hopping integrals as 
t = W/(4y/z), where W refers to the band width. In this context, an important 
research effort has been dedicated to the half-filled Hubbard model. This concerns 
mainly the metal-insulator transition within the paramagnetic phase and also the 
properties of the antiferromagnetic phase which in the absence of frustrations is the 
most stable solution at low temperatures (bipartite lattice)cl Therefore, it would 
be of considerable interest to determine the properties of the half-filled Hubbard 
model on a Bethe lattice with finite z by using accurate numerical methods. Such 
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numerical studies could be very useful, particularly in view of the possibility of 
improving the z = oo equations by introducing 1/z corrections. 

The aim of this paper is to determine several ground-state properties of the 
half-filled Hubbard model on a z = 3 Bethe lattice as a function of the Coulomb 
interaction strength U ft. For this purpose we take advantage of a property that 
Bethe lattices share with one-dimensional (ID) chains, namely, the fact that there is 
only one path between any pair of sites in the system. This characteristic allows the 
design of a simple real-space renormalization scheme in order to apply the density- 
matrix renormalization group (DMRG) methocfl. Density-matrix renormalization 
is a very powerful technique which was first proposed a few years ago in the context 
of ID spin systems. Since then it has been rapidly extended to become one of 
the leading numerical tools for the study of low-dimensional correlated quantum 
systems including recently two-dimensional lattices of finite sizeQ. The success and 
wide range of applications found by this approach rely on two main qualities: its 
high accuracy even for systems as large as a few hundreds of sites, which allows safe 
extrapolations to the thermodynamic limit, and its flexibility concerning the model 
Hamiltonian under study (Heisenberg, t — J, Hubbard, Kondo, etc). However, 
except for the work of Otsuka on the spin-1/2 XX Z HamiltonianLj, the infinite 
DMRG calculations have always been limited to ID problems. To our knowledge, 
this is the first time that the DMRG algorithm for infinite systems is applied to a 
fermion model not having the ID topology. 

The remainder of the paper is organized as follows. In the next section the main 
details of the application of the DMRG method to the Bethe lattice are given. 
Results for the ground-state energy and several spin and charge local properties are 
presented and discussed in Sec. III. Finally, Sec. IV summarizes our conclusions. 



II. DETAILS OF THE CALCULATION 

We consider the Hubbard Hamiltonian! 

H = -t 4x Cjo + u ^ "a (!) 

(i,j),(T i 

on a Bethe lattice with coordination number z = 3 (see Fig. |l|). In the usual 
notation, c| CT (cj CT ) creates (annihilates) an electron with spin a at site i, fii a = c ia Ci a 
is the number operator, t refers to the nearest neighbor (NN) hopping integral, and 
U to the on-site Coulomb repulsion. 




FIG. 1. Illustration of the Bethe lattice with coordination number z = 3. The numbers 
label non- equivalent sites. 
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Several ground-state properties of the Hubbard model are determined using the 
DMRG methodB. This is an iterative projection technique which allows to include 
the most relevant part of the ground-state wave function on a limited number of 
many-body states. The system is partitioned into several regions in real space or 
blocks between which renormalized interactions are computed. The accuracy of 
the method is controlled by the number of states m retained for the description of 
each block. In the DMRG algorithm for infinite systems, the size of the system 
is increased at each renormalization step and the properties of the thermodynamic 
limit are determined by extrapolating the succession of finite system calculations. 




FIG. 2. (a) Renormalized Bethe lattice (z = 3, see Fig. |l|) and (b) superblock renormal- 
ization procedure (y > 3). Circles represent sites which are treated exactly and rectangles 
renormalized blocks. Notice that the total number of sites N a — 3 x 2 U — 2 increases 
exponentially with the number of renormalization iterations v. 



The renormalization procedure used for the z = 3 Bethe lattice is illustrated 
in Fig. |2|. Starting from a central site i = and the two shells of its first and 
second NN's, a new shell of NN's is added at each iteration v. The total number 
sites is thus given by N a — 3 x 2 U — 2. Notice the contrast-.with the usual DMRG 
scheme for ID systems where N a increases linearly with iJB. The central site and 
its NN's (the sites on the first shell) are treated exactly at all iterations and the 
renormalizations are performed on the two branches which start at each of the 
first-shell sites (see Fig. ^J). Notice that no renormalization is actually done until 
v = 3, i.e., the 10-site system with two sites in each block is treated exactly. In 
practice, the extrapolations of the considered properties to v — > oo converge after 
v ~ 15 iterations. At this point N a ~ 2 x 10 5 . It should be however noted 
that in a Bethe lattice the surface to volume ratio does not vanish with increasing 
number of shells L. Taking into account that the number of sites in the Zth shell 
is N s = z(z — 1)' _1 (I > 1), one finds for z — 3 that 1/2 of the sites belong to 
the outermost shell, 1/4 to the first shell below the surface, 1/8 to the second shell 
below the surface and so forth. Consequently, global properties of the system such 
the average ground-state energy per site E s are dominated by the outermost shells. 
Bulk properties corresponding to the translational invariant situation have to be 
calculated locally. For instance, the ground-state energy per site E is determined 
from E — U(no^fiQi) + (zt/2)J2 a (^ia^Oa + Co a ci a ) by using the density matrix 
reduced to the central sites i = and i = 1 (see Fig. [l]). In Sec. [II it will be shown 
that the present DMRG algorithm yields accurate results for both global properties 
including the surface and local bulk-like properties. 

For the calculations we take advantage of a theorem by Lieb which states that 
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for U > the ground-state spin S of the half-filled Hubbard Hamiltonian on a 
bipartite lattice is S = \Na — -/Vb|/2, where Na and Nb are theriiumber of sites 
belonging to the two sublattices A and B (N a = Na + Nb evenp. In the Bethe 
lattice considered in this paper, even and odd shells constitute the two sublattices 
A and B (NN hopping). At the renormalization iteration v the ground-state spin 
is S = 2 ly ~ 1 so that the calculations can be performed in the subspace of maximal 
S z = S. This exact result provides an important simplification which appears to 
be crucial in order to achieve reliable results with present computer facilities, since 
for S z — S the Hilbert space is much smaller than for S z = 0. Thus, the number 
of states to kept in the renormalized blocks can be reduced significantly without 
loss of accuracy. For example for U — 0, we have performed DMRG calculations 
by using S z = and S z = S, keeping in both cases the same number of states 
to = 20 in the renormalized blocks. After the first few iterations one finds important 
differences in the ground-state energy per site, the S z — results being 2.6 x 10~ 2 t 
higher than those corresponding to S z = S. Moreover, the sum P m of the retained 
eigenvalues of the density matrix — a good criterion to estimate the quality of a 
DMRG calculationcl — follows the same trend. Indeed, for S z = S, 1 — P m is always 
smaller than 7 x 10~ 4 while for S x = 0, 1 — P m can be as large as 2 x 10 -2 , a value 
which in practice is too large for obtaining accurate results. In addition it should be 
noted that the uncorrelated limit is the most difficult case-in DMRG calculations. 
l~P rn is in fact always smaller for finite U/t than for U — Ocil. The results presented 
in this paper were obtained by keeping m = 20 states in each renormalized block. 
This would correspond to 4 4 to 3 = 2, 048, 000 possible configurations. However, only 
between 100,000-130,000 of them belong to the targeted S z = S subspace. Still, 
the dimension of the superblock density-matrix, 16m 2 = 6400, is quite important, 
which renders the computations very demanding. 



III. RESULTS AND DISCUSSION 

In Fig. H results are given for the ground-state energy per site as a function of 
U/t. 



Fig. 3 Lepetit et al. 




U/(U+4t) 



4 



FIG. 3. Ground-state energy per site of the half-filled Hubbard model on a z = 3 Bethe 
lattice as a function of the Coulomb repulsion U/t. The solid curve is calculated locally 
at the central sites from E — f7(rioffto|) + ( 2 ^/2) X^^io-^ 00 " + ^Oct^Io-) (see Fig. [j]). The 
dashed curve is the average energy E s of the complete system including the surface. Crosses 
(plus signs) correspond to even (odd) renormalization iterations. The corresponding exact 
results for (7 = are given by the open circles. 

Since in the Bethe lattice the weight of the atoms close to the surface does not 
vanish with increasing number of shells, it necessary to discern between global prop- 
erties which include surface contributions and local properties calculated close to the 
central site i — 0. The global energy E s (dashed curve) is derived by extrapolating 
the ground-state energy per site for L — » oo. For U = the DMRG calcula- 
tions yield E s = — 1.10268i, while the exact results obtained by diagonalizing the 
finite- L tight-binding matrix and extrapolating to L — > oo is E^ x = —1.103064 (see 
Appendix). The agreement seems quite remarkable since E s is given by the contri- 
butions of the atoms of the outermost shells — 1/2 of the sites belong to the surface, 
1/4 to the layer below the surface, etc. — which are renormalized already from the 
very first iterations. Similarly good results are obtained for the local ground-state 
energy E (solid curve) which is obtained from the density matrices at the central 
sites extrapolated for L — > oo. For U — 0, we obtain E = —1.52474, while the inte- 
gral of the single-particle density of states of the Bethe lattice is E ex = —1.52554. 
The renormalized blocks thus provide a proper embedding of the central sites. It 
is worth noting that the accuracy of these results, derived by keeping only m = 20 
states per block and setting S z — S at each iteration, is comparable to the accuracy 
of calculations with rn = 100-150 in the ID Hubbard model (S z = 0). It goes 
without saying that m — 100 calculations on a Bethe lattice are hardly feasible 
with present computer facilities. Since the uncorrelated limit is-the most difficult 
(less precise) case in DMRG calculations on the Hubbard modeH, we may expect 
that the results for finite U/t are at least as accurate. This is also confirmed by 
the sum P m of the retained eigenvalues of the block density-matrix which increases 
with increasing U/t. 

E and E s increase monotonously with U/t and vanish as expected in the Heisen- 
berg limit. Their U/t dependence are very similar. In fact, allowing for a rescaling 
of the energies at U = 0, E/EiJJ = 0) and_B s / 'E S (U = 0) are close to the corre- 
sponding results for the ID Hubbard chainEl. Quantitatively, E s is always higher 
than E due to surface boundary effects. Surface sites at the outermost shell have a 
smaller local coordination number z — 1. Therefore the effective local band width 
and the binding energy of surface sites are smaller than in the bulk. 

Several local properties have been calculated around the central site i = in 
order to analyze the behavior of the ground-state in the bulk limit. Results for the 
average local magnetic moment (5 Z (0)) are shown in Fig. ||. 
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Fig. 4 Lepetit et al. 
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FIG. 4. Average magnetization (^(O)) for the half-filled Hubbard model on a Bethe 
lattice (z = 3) calculated at the central site i = 0. Circles indicate tight-binding exact 
results (U = 0) or perturbation theory estimates (U = oo). The dot is the DMRG result 
for the Heisenberg model on the z = 3 Bethe lattice. 

(<Sz (0)) increases monotonously with U/t reaching 0.35 in the Heisenberg limit. 
This reflects the formation of a spin-density wave (SDW) in which the average 
spins on the two sublattices A and B point in opposite directions. Indeed, for 
NN sites (i = and i = 1) we obtain (S z (l)) ~ -(4(0)) and (4(1')) - (4(1)) 
(1(^(0) + 4(1)) | < 10~ 4 ). Notice that the calculated (4(0)) docs not vanish in 
the uncorrelated limit as it should. From the tight-binding solution for the Bethe 
lattice with a finite number of shells L one obtains (S z (0)} = 1/(3L + 2) for L even 
and (4(0)) = for L odd (see the Appendix). Thus, in finite-size systems with 
even L the local magnetization does not vanish even for U = (S z = S). Our 
DMRG calculations follow precisely the exact results during the first renormaliza- 
tion iterations but for large v a slight increase of (4(0)) is observed for L odd which 
yields a small non-zero value in the extrapolation to L — > oo ((S Z (U = 0)) = 0.025). 
We expect that this inaccuracy will be remedied by increasing the number of states 
kept in each renormalized block. In any case, since the DMRG method has better 
convergence properties with increasing U/t, the precision of the results improves 
rapidly for finite values of the interaction. In fact, in the strongly correlated limit 
the calculated (5^(0)) = 0.35 is in very good agreement with perturbative estima- 
tions and with DMRG calculations using the Heisenberg Hamiltonian. 

In order to analyze the strongly interacting Heisenberg limit we approximate the 
ground-state wave function by including first-order perturbations to the Neel state 
4>o ■ 

= </>o - J2( S t S ^ + S 7 S ^o ■ (2) 

Notice that the coefficient of the first-order correction (spin-flip states) is inde- 
pendent of the exchange constant J = 4t 2 /U, since the off-diagonal matrix ele- 
ments (J/2) and the energy differences [2(z — 1) J/2] are both proportional to J 
in the Heisenberg model. The average of local operators O [e.g., O — 4(0) or 
O = 4(0)4(1)] are obtained from 
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(3) 



where p refers to the reduced density matrix. For example, for a single site Eq. (g) 
yields p(T, T) = 1, p(i,i) = 3 /16 and p(t I) = p(lA) = 0. In this way one 
obtains (S z (i)) = 13/38 = 0.342, which compares very well with the DMRG result 
(S z (i)) = 0.348 for U/t = 128. We conclude that the reduction of (S z (i)) with 
respect to the saturation value 1/2 at very large U/t is the result of quantum 
spin fluctuations involving the first NN shell of atom i. In addition, we have also 
performed DMRG calculations for the Heisenberg model on the z = 3 Bethe lattice 
which confirm the Hubbard results at U/t > 1 [(S z (i)) = 0.347, see Fig. Q. 



Fig.5 Lepetit et al. 
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FIG. 5. Bethe-lattice results for the average spin fluctuations (S\) — (S z ) 2 as a function 
U/t, where S z = 5 Z (0) (single site, solid curve), S z = S z (0) + 5,(1) (NN spins, dashed 
curve) and S z = S z (l) + S z (l') (NNN spins, dotted curve). See Fig. oj. Circles indicate 
tight-binding exact results (17 = 0) or perturbation theory estimates (U = oo). The dots 
are DMRG results for the Heisenberg model on the z = 3 Bethe lattice. 

Fig. 6 Lepetit et al. 
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FIG. 6. Results for {^(O) 2 } of the half-filled Hubbard model as a function U/t. 
The solid curve corresponds to the z — 3 Bethe lattice and the dashed curve to the 
one- dimensional chain (z = 2). 

In Fig. ^ results are given for the fluctuation of the local magnetic moments at 
a single site and at pairs of NN and NNN sites. In all cases (S 2 ) — (S z ) 2 presents 
a maximum for U/t = 2.4-2.9. This behavior is a consequence of the interplay 
between the (S 2 ) increase which is due to the reduction of the weight of doubly 
occupied and empty sites by correlations, and the formation of permanent mag- 
netic moments (S z ), which reduces the fluctuation of the spin moments around 
their average (see Fig. ||). As shown in Fig. ||, (S 2 ) in the Bethe lattice increases 
monotonously with U ft very much like in the ID Hubbard chain. The main quali- 
tative difference between Bethe-lattice and ID results for (S 2 ) — (S z ) 2 comes from 
(S z ) which is zero for the ID case. Notice that the DMRG calculations are in 
good quantitative agreement with the tight-binding analytic results (open circles, 
U = 0), with independent DMRG calculations for the spin-1/2 Heisenberg model 
(dots, U — oo) as well as with the first-order-perturbation estimation from Eqs. @ 
and (||) (open circles, U = oo). 



Fig. 7 Lepetit et al. 
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FIG. 7. Spin correlation functions (S z (i)S z (j)) - (SM))(S z (j)) between NN sites 
(dashed curve) and NNN sites (dotted curve). As in Fig. g, circles indicate tight-binding 
exact results (U = 0) or perturbation theory estimates (U = oo) and dots are DMRG 
results for the Heisenberg model. 

The spin correlation functions (S z (i)S z (j))- (S z (i)) (S z (j)) between NN and NNN 
sites i and j are given in Fig. |7| NN spins show antiferromagnetic correlations 
which tend to become stronger as U jt increases. A shallow maximum is observed 
approximately at the same value of U/t for which (S z ) — (S z ) 2 is maximal. No- 
tice that a considerable part of the antiferromagnetic correlations between NN is 
already present for U = 0. In contrast, the spin correlations between NNN's are 
much weaker. In this case, parallel alignment is slightly favored as U/t increases 
[(S Z {1)S Z (V)) - (S z (l))(S z (V)) > 0]. These trends are consistent with the sign 
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alternations found in (S z (i)} for i belonging to different sublattices (SDW). As ex- 
pected, good agreement is obtained between the Hubbard results for large U/t and 
independent DMRG calculations for the the Heisenberg model. The first-order esti- 
mates [Eq. (0)] are somewhat less accurate in this case but still remain qualitatively 
correct. 

The average density (n(i)) at the central sites i = and i = 1 is very close to 
1 independently of U/t, which confirms the expected absence of a charge-density 
wave for U > (\(n^ +fHl) ~ 1| < 10 -4 ). The fluctuations of the density at a single 
site and at pairs of NN and NNN sites are given in Fig. |s|. 




FIG. 8. Bethe-lattice results for the density fluctuations (n 2 ) — (n) 2 as a function U/t, 
where h = n(0) (single site, solid curve), n — n(0) + n(l) (NN sites, dashed curve) and 
S z = n(l) + n(l') (NNN sites, dotted curve). n{i) — + hii, see Fig. ^. 

In all cases we observe a monotonic crossover from the uncorrelated regime 
((n 2 ) — (h) 2 maximal) to the strongly correlated, localized regime where charge 
fluctuations are suppressed. The U/t dependence is quite similar to what is ob- 
tained for ID Hubbard model. The density fluctuation at a pair of NNN sites 
(dotted curve) is approximately twice the single-site result (solid curve) which 
indicates, as in the case of the spin degrees of freedom, that density correla- 
tions between NNN's are very weak [(n(l)n(l')) — (n(l))(n(l')) ~ 0]. In con- 
trast, for NN sites charge fluctuations are significantly smaller. One observes 
that ([n(0) + n{\)] 2 ) - (n(0) + n{\)) 2 ~ (3/2)[(n(0) 2 ) - (n(0)) 2 ] or equivalently 
(n(O)n(l)) - (h({)))(h(l)) ~ [(n(0) 2 ) - (n(0)) 2 }/2. Notice that these relations hold 
approximately for all values of U/t, even for U = 0. The ratio between single- 
site charge fluctuations and fluctuations on a pair of NN's is not much affected by 
changes in the Coulomb repulsion strength and therefore seems to result mainly 
from the geometrical proximity of NN sites. 
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IV. CONCLUSION 



Several ground-state properties of the half-filled Hubbard model have been de- 
termined on the Bethe lattice with coordination z = 3 by using a density-matrix 
renormalization group (DMRG) algorithm for open infinite systems. Though the 
lattice is not one dimensional (ID), the existence of a unique path between any pair 
of sites allows to formulate a simple renormalization procedure. In contrast to pre- 
vious density-matrix renormalizations studies of Hubbard-like models on ID chains 
or ladders where the number of sites N a increases linearly with the number of itera- 
tions v, in the present approach N a increases exponentially with v (N a — 3x2^ — 2). 
This is a consequence of the fact that 2 blocks are renormalized into a single one 
at each iteration. Despite the very rapid increase of N a , the DMRG method pro- 
vides accurate results over the whole range of U/t already by keeping few states 
per block (m = 20). This is achieved by working in the subspace of maximal spin 
projection S ZJ j= S, where the ground-state spin S = 2 U ~ 1 is derived from a the- 
orem by Liettj. For example, in the limit of U = the calculated ground-state 
energy per site differs by only 3 x 10 _4 t from the exact tight-binding result. It is 
remarkable that this level of precision concerns not only local properties calculated 
at the unrenormalized central sites, but also global properties which are dominated 
by the renormalized sites of the outermost shells. From a general point of view, the 
present study encourages renormalizations of more than one block into a superblock 
in future DMRG algorithms. 

The main results for the z — 3 Bethe lattice may be summarized as follows. The 
calculated local magnetic moments (S z (i)} increase monotonically with increasing 
Coulomb repulsion U/t forming an antiferromagnetic spin-density-wave state which 
matches the two sublattices of the bipartite Bethe lattice. The maximum (S z (i)) 
found in the Heisenberg limit ((S z (i)) — 0.35) is reduced with respect to the sat- 
uration value (S z (i)) = 1/2 as a result of exchange flips between the spin at site i 
and its NN's that point in opposite directions. The fluctuations of the local spins 
(S z (i) 2 ) — (S z (i)) 2 show a maximum as a function of U/t for U/t — 2.4-2.9. For 
small U/t the usual increase of (S z (i) 2 ) due to the reduction of double occupa- 
tions dominates, while for large U/t the formation of large permanent moments 
(S z (i)) blocks local spin fluctuations. NN sites show important antiferromagnetic 
spin correlations which increase with increasing Coulomb repulsion. However, NNN 
sites are very weakly correlated over the whole range of U/t. The AF correlations 
are very short ranged in contrast to the static picture of a SDW. This reflects the 
importance of quantum fluctuations in the z = 3 Bethe lattice as in ID systems. 

Taking into account that the Bethe lattice is one of the standard models for the 
studying the limit of infinite dimensions, it would be very interesting to compare 
the present results for z = 3 with the outcome of the d = oo equations in order to 
quantify the importance of X/d corrections. Moreover, further DMRG investigations 
including NNN hoppings should be very valuable for a better understanding of the 
metal-insulator transition in the presence of frustrations, particularly since the in- 
formation derived from the DMRG method is complementary to finite-temperature 
quantum Monte-Carlo calculations. 
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APPENDIX 



The aim of this section is to outline the block diagonalization of the tight-binding 
matrix of a finite Bethe lattice formed by a central site i — and L successive 
nearest neighbor (NN) shells. The number of sites at shell I (1 < I < L) is given by 
N s (l) — z[z — 1) , where z refers to the coordination number. The notation used 
for labelling the lattice sites is illustrated in Fig. ^|for z = 3. 




FIG. 9. Labelling of Bethe-lattice sites used in this appendix. 



A site belonging to the shell I is denoted by the set of I numbers (ii, «2> • • • i;) 
which define the path to follow in order to reach the desired site starting from the 
central site i — 0. The tight-binding matrix Hq of the z = 3 Bethe lattice with 
NN hoppings t can be block diagonalized by using that the finite Bethe lattice is 
invariant after the transposition of any pair of branches that connect a site of shell 
I — 1 with its 2 NN's of shell I. For the outermost shell (7 = L) the symmetry 
adapted single-particle states are 



|ii,i2) ■ ■ -«l-i,+) = {\h,i2, ■ ■ -ii-i, 1) + |*i, k, ■ • .i£-i,2)) /V2 
|zi,i 2 , . . . IL-1, -) = *2, - - - iL-i, 1) - \h,i2, ■ ■ ■ ii-i,2)) /v2 . 



(4) 
(5) 



For the other shells (I < L) on proceeds recursively in the same way building 
symmetric and antisymmetric linear combinations. For the first shell the 3- fold 
symmetry around the central site i = is applied. 

In the new basis Hq splits into a (L + 1) X (L + 1) matrix block of the form 



y/3t 
Vit V2t 
V%t y/2t 



V2t 

Vzt o 



(6) 



which involves only purely even states including the central site, and in smaller / x I 
matrices Bi with 1 < I < L of the form 
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B, = 



V2t 

V2t o Vzt 

y/2t y/2t 



V2t 
V2i 



(7) 



£?l appears twice in iJ an d each of the other Bi appears 3 x 2 L ~' _1 instances in 
the total tight-binding matrix. 

The eigenvalues of Bi are (3k — —2^/2tcos ^ with k e [1, Z] and those of A arc 

afe = —2^21^1 cosO^ where the 9k are the roots of 2sin(L + 2)# = sini^. The 
later equation is solved numerically and the tight-binding energy per site E s of the 
finite Bethe lattice is determined for any L. The extrapolated value for L — > oo 
is E s = — 1.10306i. Notice that this result differs from the integral of the local 
density of states at the central site (E — — 1.5255i) since E S (L — > oo) is dominated 
by surface contributions. 
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